Ground-state phases of ultracold bosons with Rashba-Dresselhaus spin-orbit couphng 
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We study ultracold bosons in three dimensions with an anisotropic Rashba-Dresselhaus spin- 
orbit coupling. We first carry out the exact summation of ladder diagrams for the two-boson t 
matrix at zero energy. Then, with the t matrix as the effective interaction, we find the ground-state 
phase diagrams of bosons in mean field, as a function of the spin-orbit coupling, the anisotropy, 
and the scattering lengths between particles in the same and in different pseudospin states. The 
resulting phase diagrams have much richer structures than those obtained using mean-field couplings, 
exhibiting three different phases: a plane wave condensate, a striped condensate, and an unstable 
phase. The differences between the present approach using the t matrix compared to using mean- 
field couplings is significant for large scattering lengths, large spin-orbit-coupling strength, or small 
anisotropy. 



I. INTRODUCTION 

Spin-orbit coupling plays a crucial role in a variety 
of physical systems ranging from atoms and nuclei to 
topological insulators [ill and spintronics • Recently, 
the prospect of realizing spin-orbit coupling in ultracold 
atomic systems has led to increased interest not only in 
fermions but also bosons for this purpose Q. Currently 
proposed and realized schemes to produce spin-orbit cou- 
pling in ultracold atomic systems create effective non- 
Abelian gauge fields for atoms which take the form 
of Rashba-Dresselhaus smn-orbit couplings, familiar in 
semiconductor physics [lol . ITl| . A system of spin-orbit 
cou;^ed ultracold bosons was recently realized by Lin et 

The problem of ultracold bosons with Rashba- 
Dresselhaus spin-orbit coupling has been considered 
within mean- field, where the interparticle interactions are 
assumed to be independent of momenta [l^ - [T7| . These 
studies predicted that the ground states can be either 
a "plane-wave" Bose-Einstein condensate (BEC) of par- 
ticles in a single momentum state, or a "striped" BEC 
involving a coherent superposition of two different mo- 
menta. However, recent studies indicate that, beyond 
mean field, the effective interaction has a qualitatively 
different structure, resulting, for example, in the absence 
of interaction between particles scattering in the zero 
total-momentum channel [18, 19] and the prediction that 
a BEC involving only a single momentum of particles is 
not favorable in an isotropic Rashba field [l^ . 

In this study we delineate the ground state phase 
diagram of bosons taking the effective interaction de- 
scribed in terms of the full t matrix rather than mean- 
field coupling. We consider bosons with an isotropic 
or anisotropic Rashba-Dresselhaus spin-orbit coupling, 

K,{axPx + V^yPy): where n is the spin-orbit coupling 
strength, p is the particle momentum, and 77 deter- 
mines the x-y anisotropy. In our previous paper p^ . 
we showed that the effective interaction can be renor- 
malized in terms of physical scattering lengths and does 
not depend on the ultraviolet cutoffs. Here, we extend 



this calculation, and carry out the exact summation of 
ladder diagrams for the t matrix of bosons scattering at 
zero energy in three dimensions. Then we reduce the 
Hamiltonian, with the t matrix as the effective interac- 
tion, to a Nozieres model [20| . from which we determine 
the ground-state phases as functions of the anisotropies 
in the naaa-^OLab plane, where aaa and aab are the scat- 
tering lengths between the same and different pseudospin 
states, respectively. We show that there are generally 
three phases: a plane wave BEC, a striped BEC with two 
different momenta, and an unstable phase where the ef- 
fective interaction is attractive [25]. The phase dia grams 
(Figsl2]and[3]below) are substantially different and richer 
in structure than those predicted using mean-field cou- 
plings with the bare couplings replaced by AirlPaij /m, 
especially when nuaa and KUab are large, or 77 is close to 
unity. In the isotropic limit rj — 1, the plane wave BEC 
does not appear (Fig.|31). We describe how the phase di- 
agram evolves continuously with anisotropy, from rj = 
where spin-orbit coupling is present only in the x direc- 
tion, to an isotropic spin-orbit coupling r] = 1. We find 
that, in the vicinity of rj = 1, the phase diagrams are log- 
arithmically sensitive to small changes in 77, a structure 
we analyze by expanding in the anisotropy about 77 = 1. 
We also find that a BEC with negative scattering lengths 
can be stabilized in the presence of spin-orbit coupling, 
when the scattering lengths are sufficiently large in mag- 
nitude (Fig. 

Our analysis can be directly compared with proposed 
experimental schemes when the scattering lengths 
are all equal; a good approximation for the three F = I 
hyperfine states of *^Rb. We discuss in this situation 
how the phases of bosons in the ground state evolve with 
varying anisotropy and scattering length (Fig.[S]). In the 
following we take h = 1. 



II. HAMILTONIAN 

We consider bosons in three dimensions with a Rashba- 
Dresselhaus spin-orbit coupling, described by the Hamil- 
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tonian 



the rj dependence: 



^ = E (4 ^1.) 
p 



(5aa'^P4 ^P2 '^Pl 



Pl+P2=P3+P4 
+5bb^P4 ^P3 bpi + '^dab a],, ^j,3 ^P2 flpi ) ■ ( 1 ) 



As in our previous paper [19|, m is the atomic mass, V 
is the volume of the system, Cp annihilates an atom the 
pseudospin state a with momentum p, and bp annihilates 
an atom in the pseudospin state h with momentum p; the 
(Jx and Gy are the usual Pauli matrices between the in- 
ternal states, and / is the two-by-two identity matrix. 
We take the coupling k to be positive. The Qaa, gbb, and 
gab are the bare s-wave couplings between a-a, b-b, and 
a-b particles. When the system is isotropic, rj = I, the 
Hamiltonian completely reduces to that considered pre- 
viously [19] . In practice the effective Hamiltonian, in the 
basis in which the coupling has the Rashba-Dresselhaus 
form, can contain terms that do not conserve the individ- 
ual number of particles of each species (a- like and 6- like) ; 
we ignore such terms at this point, and return to this is- 
sue below. 

To diagonalize the single-particle part of the Hamilto- 
nian we introduce operators ap and /3p defined by 



Pl+P2=P3+P4 
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+ v 
+ v. 
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(5) 

il,02;03,04 



('^P4 '^P3 '^P2 '^Pl 

(/3p4/3p3"p2api +"p4"P3/^P2^Pi) /2 

("P4^P3/^P2/3P1 +/3p4"p3"P2"Pl) 
("P4"P3^P2"P1 +^P4/^P3"P2^Pl) 
'^P4/^P3/^P2'^Pl 



(5) 



where 4>i is the angle of iPi,x,VPi,y) in the x-y plane, 
dependent on rj, and the V*^*^'s are 



(1) 

1,02;<P3,<P4 

(2) 

<t>i:4'2;4>3,4>i 
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1,02;03,04 
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t>l,4>2;<p3,'l>4. 
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with 



(2) 



where (/) is the angle of {Pxt'HPv) the x-y plane. Note 
that, for a given momentum, (j) depends on 77. In terms 
of the as and /3s, the Hamiltonian becomes 



H = Y, (e-(p)«^,«P + e+(p)/5^/3p) + ^i: 



The single-particle spectrum has two branches 



2m 



(3) 



(4) 



i{4>l+<t>2-4>3~'t>4) 



)/4. 



(7) 



III. THE f-MATRIX 



We now turn to calculating the scattering t matrix, 
rj^^(p, p'; q), for two bosons in the ground states of the 
lower (a) branch, with incoming momenta q/2 -I- p and 
q/2 — p, and outgoing momenta q/2 -I- p' and q/2 — p'. 
Note that, for particles in the single particle ground state, 
the z component of the momentum is zero. As we showed 
earlier, the ultraviolet divergences in the t matrix can be 
renormalized in terms of cutoff-independent low-energy 
parameters p^] . and, in the isotropic case, the t matrix 
depends only on the spin-orbit coupling strength, the 
scattering lengths, and the incident and outgoing mo- 
menta. The exact summation of ladder diagrams for the 
t matrix, for both isotropic and anisotropic spin-orbit 
couplings, is then (see Appendix A for the derivation): 



the single particle ground state is given by the lower 
branch e-(p), which has degenerate states on the cir- 



cle y + = K and Pz = when rj ^ 1, and twofold 
degeneracy for p = (±k, 0, 0) when < 77 < 1. 

The interaction part Hint is as in Ref. [l^, except for 



no(p,p';q) 



TT / -J- f.*'f'2 
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where 
M = 

fm + 

hm/2) 



/ii(q/2) 

f(i) + ^ 

h2iq/2) 



fii) 



/i2(q/2) 
/i^(q/2) 



(9) 



with q = q/n and g = g/K. The angles (jji, (1)2, 4'3i ^nd 
(j>4 are those of the vectors q/2 + p, q/2 — p, q/2 — p', and 
q/2 + p' in the x-y plane with a factor of 77 multiplying 
the y components. The dimensionless functions /(q/2), 
g{q/2), /ii(q/2), and /i2(q/2) are, 
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where (^5 and 4>q are the angles of q/2 — k and q/2 + k in 
the x-y plane with y components multiplied by rj. Chang- 
ing the angle of q in the x-y plane only changes the overall 
phases of /ii(q/2) and /i2(q/2). These four functions are 
everywhere finite except for the logarithmic divergence 
of f{q/2) at g = 0. When the scattering lengths are 
small, the diagonal elements of M are dominant and we 
may ignore the off-diagonal elements; the t matrix thus 
obtained does not have terms containing products of dif- 
ferent scattering lengths, which is the _approximate result 
obtained for the bosonic t matrix in [19i . 



IV. GROUND STATE PHASES 

We now determine the many-body ground state via 
mean-field theory using the t matrix derived above as 
the effective interactions, an approximation valid as long 
as the nafj are all ^ 1, where n is the particle density. 
In mean field, we assume that all particles are in the 
single-particle ground states (k, 0,0) or (— k, 0,0), and 



thus ignore possible occupation of excited states as a 
consequence of the interaction. In this case the system 
is described essentially by the Nozieres model [l^]. The 
issues of going beyond mean field, (e.g., via Bogoliubov 
theory) as well as including possible effects of the conden- 
sate on the effective interaction, are beyond the scope of 
this paper and are left for the future. For < 77 < 1, 
we take the particles to be either at p = (k, 0, 0) or 
(— K,0, 0); the relevant interactions are those between 
particles of either the same momentum or opposite mo- 
menta. We denote the interaction with same momentum 
by To = r""(0, 0; ±2K) and that with opposite momenta 
by r^r = r""(±K, =fK, 0), where we use the abbreviated 
notation K = (k, 0, 0). 

The relevant terms in the interaction are then equiva- 
lent to the Nozieres model ^] 

^int ~ ^roiV(7v + 1) + -i(2r, - ro)7V.7Vo, (n) 



"(«:,0,0)"(«^.0,0) 



and A^TT = 01] 



K,o,o)"(-«^o,o)- 



where A^o = "(k,o,o)"('^.'J,u) f^^^" ^'tt = 
The total number of particles, N ~ Ni^ + N^^, is fixed. 
For Tq < 2r^, the ground state is a single BEC with 
either all the particles in (k,0, 0) or (— k, 0,0), while, for 
Fq > 2F7r, the condensate is nominally fragmented with 
half of the atoms forming a BEC in one state and the 
other half forrning a BEC in the other state. However, as 
shown in Ref. [2J| , such a fragmented state is expected to 
be unstable against formation of a coherent condensate 
with a condensate wave function that is a coherent super- 
position of the two momenta. Following the conventions 
of Refs. [13, [l^, we call the single BEC phase "plane 
wave," and the BEC phase with two different momenta 
"striped." The difference of the present calculation from 
earlier studies with mean-field couplings [3, [iBl , is that 
here the bare couplings, Vq^o-q^o 

and V^]l 

;0,7r: replaccd 

by Fq and F^, respectively. 

While there is no difficulty in deriving the phase di- 
agrams for general scattering lengths, we assume here 
for simplicity that the intraspecies scattering lengths are 
equal, Gaa = abb- Then, 

2tt 

Fo = X 

rriK 

l/Kagg + l/naab + 2/(1) - g(l) + hijl) - 4/12(1) 

{1/Kaab + /(I) - g{l)){l/Kaaa + /(I) + hi{l)) - 4/7,2(1)2 

r I (12) 

where /ii(0) = /ii(q = (0,0,0)), hi{l) = 
= (1,0,0)), etc. The quantities 

/(0),/7i(0),/(l),g(l),/ii(l), and /i2(l), which depend 
on 77, can be calculated numerically. The interaction 
between different momenta F^ is independent of aab, 
and is a monotonically increasing nonnegative func- 
tion of HMaa, cqual to at naaa = and reaching 
27r/ [7?7.k(/(0) — /ii(0))] at naaa = 00. The dependence 
of Fq on Kflaa and naab is more complicated. We plot Fq 



4 



and r^r, both scaled by 27r/(mK), for ij = 0.5 in Fig. [T] 
Now we discuss the ground-state phases from rj — to 1. 




1.0 



(a) 




(b) 



FIG. 1: (Color online) (a) Fq as a function of naaa and KUab, 
and (b) F^ as a function of KUaa, both scaled by 27r/(m«:), 
for rj — 0.5. The vertical plane in the middle of panel (a) 
indicates the resonance where, from left to right, Fq diverges 
to positive infinity and comes back from negative infinity. 



When 1] = 0, the effective interactions are relatively 
simple. It can be shown that /(O) = /ii(0) for = 0; 
hence F^ — 2-Kaaa/fTii and the effective interaction in 
the q = channel does not depend on the spin-orbit 
coupling strength k. In the q/2 — (k, 0, 0) channel, 
/(I) = -1,5(1) = 0,/ii(l) = 0, and h^il) = 1/2, so. 



Fn = 



27r Kan 



KUab - "iKUaaKaab 



niK 



1 f^^aa I^O^ab 



for T] = 0. The effective interaction at small nuaa and 
KUab is positive, and diverges when nuaa + nuab ap- 
proaches unity. As one crosses the line Kaaa + naab ~ 1 , 
Fq starts at negative infinity and remains negative un- 
til = Anaaa^aab, after which Fq is posi- 
tive. When Fq is negative, we expect the BEC in bulk to 
be unstable against collapse, as in ordinary BECs with 
negative scattering length in the absence of spin-orbit 
couplings. We call the phase with an attractive interac- 
tion "unstable." The three possible ground-state phases. 



plane wave, striped, and unstable, are determined by the 
sign of Fq and the interplay between Fg and F^r. 

As rj increases from 0, the basic structure of Fq does 
not change; Fq remains positive at small naaa and ^^aab^ 
and as these variables increase, Fq again diverges at a 
line in the naaa-no-ab plane, beyond which it is negative 
up to a second line, after which Fq is positive. Since the 
denominator of Eq. (1121) for Fq is quadratic in l/(Ka), 
it has in fact two zeros, the one for positive scattering 
lengths, as shown, and a second for negative scattering 
lengths, which is discussed at the end of this section. The 
structure for positive scattering lengths is illustrated in 
Fig.[Tl for T] = 0.5. 
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FIG. 2: (Color online) Ground state phase diagrams in the 
naaa - K,aab plane for anisotropies = 0, 0.25, 0.5, and 0.75. 
The regions P are the plane-wave phase with a BEC of a 
single momentum. The regions S are the striped phase with 
a BEC of a coherent superposition of two different momenta. 
The phase in the regions U are unstable, with the effective 
interaction Fq negative. Along the line between S and U, Fq 
diverges, and along the line between U and P, Fq vanishes. 
The intersection of these two lines is a critical point. The 
dashed lines indicate the phase diagram derived using mean- 
field coupling, in which the plane is separated into an upper 
striped region and a lower plane wave region. 



The ground-state phase diagrams for various rj are 
plotted in Fig. [2] In the panels, the plane-wave phase 
is labeled "P," the striped phase "S," and the unstable 
phase "U." The plane- wave phase occurs when Fq < 2F^, 
the striped phase when Fq > 2F^, and the unstable phase 
when Fq < 0. Note the overall tendency of the phase di- 
agrams as rj increases; the upper striped region detaches 
from the resonant critical point, where the resonant line 
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(between S and U) and the line with Fq = (between 
U and P) touch, and the region is pushed upward as 77 
is increased. Meanwhile, the shapes of the resonant line 
and the boundaries of plane- wave regions change but, 
with the exception of the upper striped region, the overall 
topology does not change. The dashed lines 
the figures are the phase-separation lines obtained earlier 
using mean-field couplings ATraaa/'m and 47raaf,/m; 
there the striped phase is preferred above the dashed lines 
and the plane-wave phase is preferred below the dashed 
lines. Use of mean-field couplings is accurate for small 
KGaa and KGab but, as these variables increase, the de- 
viation from the mean- field- coupling prediction becomes 
significant and the phase diagrams exhibit qualitatively 
new and rich structures. 

This overall tendency continues to around rj ^ 0.99. 
With further increase of r/ toward isotropy, ry = 1, we 
start to observe qualitatively new behavior of the phase 
diagrams. The phase diagrams close to 77 = 1 are plotted 
in Fig. [21 As one sees, the striped region comes back from 
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/ii(l), and h,i{V)i to their values at 5 = and approx- 
imating /(O) by |ln(5|/4, corresponds to fixing Fp and 
varying the slope of Ft^. In the isotropic limit (5 = 0, 
F„ = and thus a plane-wave region is not allowed [cf. 
Eq. (fTTj) ]. With small anisotropy, F^ can be positive, and 
small plane- wave regions appear. 

We briefly consider tuning the scattering lengths to 
negative values. In the absence of spin-orbit couplings, 
negative scattering lengths lead to an instability in large 
systems. On the other hand, as we see from Eq. (IT2|) . tun- 
ing the inverse scattering lengths to just below does not 
immediately lead to an attractive interaction; in the pres- 
ence of the spin-orbit coupling fields, Rashba-Dresselhaus 
couplings can stabilize BECs with negative scattering 
lengths if the inverse scattering lengths are small. Even 
when Fq is negative, systems with small particle num- 
ber can be metastable in the presence of an attractive 
interactioni26]. For illustration, we plot the phase dia- 
gram for ?7 = 0.5 extended to negative scattering lengths 
in Fig. m In the regions marked "Stable," Fq > and the 
ground state is either a plane-wave or striped phase. As 
seen in the figure, when both scattering lengths aaa and 
flab are negative and large, another stable region appears 
in the phase diagram, in which the ground state is in the 
striped phase. The line between the lower-left striped 
phase and the unstable phase is a second resonant line 
along which Fq diverges. A stable region with negative 
scattering lengths generally exists for all < ry < 1; as 
increases, the stable region in the phase diagram becomes 
larger. 
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FIG. 3: (Color online) Ground-state phase diagrams for r\ 
close to unity. 
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above and touches the resonant critical point, and at the 
same time the lower plane- wave region detaches from the 
critical point. In the limit 77 = 1, the plane- wave region 
vanishes. 

The behavior around 77 ~ 1 is in fact logarithmic in 
the deviations of the anisotropy r\ from unity. We write 
5 = 1- ry2; as J ^ 0, /ii(0), /(I), 5(1), h^{\), and h^{\) 
approach finite values, but, in leading order for small (5, 
/(O) - |ln5|/4. Setting, for small 5, /ii(0), /(I), g(l). 



FIG. 4: (Color online) Ground-state phase diagram for r\ = 
0.5 extended to negative values of scattering lengths. The 
regions marked U and S are unstable and striped phases, as 
before. The region marked "Stable" is either a plane-wave 
or striped phase. Note the appearance of a stable (striped) 
phase when both scattering lengths are large and negative. 
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V. CONCLUSION 

Proposed schemes to realize Rashba-Dresselhaus spin- 
orbit couplings in ultracold atomic experiments y-|8| 
use Raman lasers to couple atoms in different hyper- 
fine states. In general as one transforms the original 
basis to one in which the coupling has the Rashba- 
Dresselhaus spin-orbit structure, the interaction Hamil- 
tonian acquires terms such as <ip^(ip^o,p2bpi which do not 
conserve the number of particles in each pseudospin state 
(a-like and 6-like). Our analysis, which did not take such 
terms into account, can be directly compared with pro- 
posed experiments when the interaction is independent 
of species (gaa = 9bb = gab), in which case the interaction 
is independent of the choice of basis. This condition is 
a good approximation for the three hyperfine states of 
^''Rb in the lowest F — 1 state. The assumption that 
gaa = gbb — gab corrcsponds to the (dashed) diagonal 
lines in Figs. [2] and El Figure O shows the phase dia- 
gram in the rj-Ka plane, where a is the assumed common 
scattering length. 
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FIG. 5: (Color online) Ground state phase diagram when 
cLaa ~ a-bb = clab — am the rj-Ka plane for (a) anisotropics less 
than 0.99 and (b) anisotropics close to unity. The horizontal 
axis of panel (b) is a logarithmic scale. 



0.9 < ?7 < 1, as drawn on a logarithmic scale in Fig. [5Kb), 
we find that the line separating the lower plane- wave and 
striped regions terminates and another line starts from 
positive infinity above which the striped phase is pre- 
ferred. This new line touches the uppermost line (below 
the upper P phase) in the figure in the 77 — > 1 limit and 
thus no plane-wave region exists at isotropic spin-orbit 
coupling [2 7]. 
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Appendix A: Exact summation of ladder diagrams 
for the t-matrix 



Here we outline the derivation of the exact t matrix ([5]) . 
Our starting point is the set of Bethe-Salpeter equations: 



rS"(p,p';q) = v«^^^. 



>3,<P4 



V^^^ F^'^l'k n'- n\ 

e_(f-k)-f6_(|+k) 



+ ,^(|_k) + e+(H+k) 
r^)^(p,p';q) = "l^0i^02;03,0^ 



-k) + 6-(f +k) 



(2^ 



V, 



(2) 



/gg(k,p^q) 
e_(f-k)-f6_(§+k) 

,(3) 



and 



'e+(f-k) + e+(§+k) ^ e+(f-k) + e_(§+k) 



rs?(p,p';q) = vitl^,^^3, 



d^k 
(2^ 



V 



(4) 



fe.0/gg(k,P^;q) 
,_(3_k)+._(|+k) 



V£0.;0.>.rgg(k,p^q) 



V 



i(5) 



,^rgg(k,p^q) 



e+(f - k) + e+(f k) + e+(f - k) + 6_(§ + k 



(Al) 



For < 77 < 0.99 the system, with increasing Ka, ex- 
periences transitions from plane-wave to striped, then to 
unstable, and finally to the plane-wave phase again, as 
seen in Fig. [SJa). Looking more closely at the region 



where r^^(p, p'; q) is the t matrix for scattering of par- 
ticles in the branches fi^v with momenta q/2 ± p to 
branches p, t with final momenta q/2 ± p' . The key 
to solving this set of equations is to construct the quan- 



7 



titles Rewriting the Bethe-Salpeter equations in terms of X, 

Y, and Z, we see that X, Y, and Z do not depend on 



1 



-'^(P,p';q) = T (r"Q(p,p';q) +r^;;^(p,p';q) their first arguments. Namely, we can write X(p, p'; q) = 

X{p';q), etc. Then X ,Y, and Z inside the integrals can 
+r^^(p,p';q)/V2 + r"^(-p,p';q)/%/2') , be moved outside, and we can solve for X, Y, and Z 

algebraically, and reconstruct FJjJ" from X , Y, and Z via 

i'(p,p';q)e-<»-+»'l = - (p, p'; q) + r|g(p, p'; q) 

-r;j,,p';q,/v^-rs,-p,p';q,/v^) ^ ^ 

and + Z(p'; q)(e^'^i + e*'^^)(e-*'^3 - e"*"^*). (A3) 

Zfp p'.q) ^ i rgg(p,p';q)-rgg(p,p';q) 

1 paaCr. r.'-«W /9 - ra"('-r. r,'-r,V /9 Introducing the free field scattering lengths by 
= iia£lI^P_^qj^v^ ^a/3l P,P,qj/V^ m/(47ra,j) = l/.g,j + mA/i2n^), where A is the high mo- 

2 (e'-^i - e^<t'2)(e--<P3 + e-''>4) ' mentum cutoff, we obtain Eq. ®. 
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